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(3 ■ ABSTRACT 

. Aims. In this work we present a revised measurement of the phase-averaged optical polarisation of the Vela pulsar (PSR B0833- 
45), for which only one value has been published so far (Wagner & Seifert 2000). 

Methods. Our measurement has been obtained through an accurate reanalysis of archival polarisation observations obtained 
with the FORSl instrument at the VLT. 
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O " Results. We have measured a phase-averaged linear polarisation degree Pl = 9.4% ± 4% and a position angle 9 — 146° ±11°, 

^ very close to the ones of the axis of symmetry of the X-ray arcs and jets detected by Ciandra and of the pulsar proper motion. 

^2 ' Conclusions. We have compared the measured phase-averaged optical polarisation with the expectations of different pulsars' mag- 
' netosphere models. We have found that all models consistently predict too large values of the phase-averaged linear polarization 

^ , with respect to the observed one. This is probably a consequence of present models' limitations which neglect the contributions 

! ' of various depolarisation effects. Interestingly, for the outer gap model we have found that, assuming synchrotron radiation for 

rS the optical emission, the observed polarisation position angle also implies an alignment between the pulsar rotational axis and 

^ . the axis of symmetry of the X-ray arcs and jets. 
_ _ I 
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1. Introduction 

Polarisation measurements of pulsars and of their syn- 
chrotron nebulae are uniquely able to provide deep in- 
sights into the highly magnetised relativistic environment 
of young rotating neutron stars. Besides the radio band, 
optical observations are primarily suited to provide such 
insights. Being the first, and the brightest {V ^ 16.5), 
optical pulsar identified so far, it comes as no surprise 
that the first optical polarisation measurements were ob- 
tained for the Crab pulsar (PSR B0531+21). Strong po- 
larisation is actually expected when the optical emission 
is produced by synchrotron radiation which, in the case 
of the Crab pulsar, is certified by its power-law spectrum 
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(see, e.g., Sollerman et al. 2000). Indeed, the polarisa- 
tion of the Crab pulsar optical emission was discovered 
(Wampler et al. 1969) soon after its optical identifica- 
tion (Cocke et al. 1969). Accurate optical phase-resolved 
polarisation measurements of the Crab pulsar (e.g. Smith 
et al. 1988; Kanbach et al. 2005) have shown the phase 
dependence of the polarisation degree, with a maximum 
of ~ 50% in the bridge between the main and the inter- 
pulse. Optical polarisation measurements are extremely 
useful in providing observational constraints to test the 
different models proposed for the magnetospheric activity 
of pulsars. Unfortunately, despite the number of pulsars 
detected in the optical band has increased significantly in 
the last ten years (see Mignani et al. 2004 and Mignani 
2005 for updated reviews), the Crab pulsar is still the only 
one for which both precise and repeated polarisation mea- 
surements have been obtained. For the other young pulsars 
(age < 10 000 years) only preliminary or uncertain phase- 
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averaged optical polarisation measurements were reported 
so far. For the second youngest 1600 years) optical pul- 
sar PSR B0540-69 {V - 22.5), considered the Crab pulsar 
"twin" because of its overall similarities in age, period 
and energetics (e.g. Caraveo et al. 2001a), polarisation 
observations with the VLT were reported by Wagner & 
Seifert (2000). However, the polarisation measurement (« 
5% with no quoted error) was certainly contaminated by 
the contribution of the compact (~ 4" diameter) surround- 
ing synchrotron nebula (Caraveo et al. 2001a). For PSR 
B1509-58, the next youngest (^ 5000 years) pulsar with 
a candidate optical identification, the value of the optical 
polarisation is also very uncertain. The original counter- 
part (V = 22) proposed by Caraveo et al. (1994) was dis- 
carded by Wagner & Seifert (2000) on the basis of the lack 
of evident optical polarisation in favour of a much fainter 
[R ^ 26) object hidden in the PSF wings of the star, 
whose detection has been now confirmed by IR observa- 
tions (Kaplan & Moon 2006). However, both the difficul- 
ties in the PSF subtraction and the faintness of the coun- 
terpart make the reported polarisation measurement (« 
10 %, also quoted with no error) only tentative. A prelim- 
inary polarisation measurement was reported by Wagner 
& Seifert (2000) also for the Vela pulsar (PSR B0833-45), 
the third brightest {V - 23.6) and youngest (- 10 000 
years) optical pulsar. Finally, for the older (~ 100 000 
years) PSR B0656-I-14 {V ~ 25) phase-resolved optical 
polarisation measurements have been obtained by Kern 
et al. (2003). Like for the Crab, also for PSR B0656-I-14 
the polarisation degree is phase-dependent, with a maxi- 
mum (~ 100 %) in the bridge between the two pulses. 
For the Vela pulsar, Wagner & Seifert (2000) pubhshed 
only the value of the linear polarisation degree. The need 
to obtain a more complete description of the polarisation 
properties of Vela pulsar by including e.g., the linear po- 
larisation position angle has prompted us to undergo a 
careful reanalysis of the original data set. The knowledge 
of the polarisation position angle is indeed an additional 
important piece of information to test different models 
on the optical emission of pulsars. This paper is organ- 
ised as follows: observations, data analysis and results are 
described in Sect. 2, while the comparison with existing 
pulsar emission models is discussed in Sect. 3. 

2. Observations and Data Analysis 

2.1. Observations 

Phase-averaged polarimetry observations of the Vela pul- 
sar were obtained at the first Unit Telescope (UTl) of 
the ESO Very Large Telescope with FORSl on April 
12 1999 as a part of the instrument commissioning. 
FORSl (FOcal Reducer Spectrograph), a multi-mode in- 
strument for imaging and (long-slit/multi-object) spec- 
troscopy equipped with polarimetric optics, is described in 
Appenzeller et al. (1998). Observations were performed in 



Table 1. Summary of the polarimetric observations of the 
Vela pulsar field taken with FORSl. 



Seq. 


Angle(°) 


Exposure (s) 


Seeing (") 


Airmass 


1 





1000 


0.81 


1.37 


2 


22.5 


1000 


0.78 


1.44 


3 


45. 


1000 


0.57 


1.53 


4 


67.5 


1000 


0.62 


1.63 


5 





380 


0.61 


1.76 


6 





380 


0.62 


1.90 



Standard Resolution (SR) mode (O'.'2/pixel) through the 
R Bessel filter (A = 6570 A; AA = 1500 A), with an aver- 
age seeing of ^ 0'.'7 and airmass varying between 1.4 and 
1.9. The instrument set-up was chosen to allow for linear 
polarisation measurements. A total of 6 frames were taken 
with the instrument position angle (hence the WoUaston 
axis) oriented at +10° North to East. One 1 000 s plus two 
380 s exposures were obtained with the fast axis of the 
half-wave retarder plate parallel to the WoUaston prims 
(rotation angle 0°) and the remaining three (1 000 s each) 
at rotation angles 22.5°, 45°, and 67.5°, respectively (see 
Table 1 for a log of the observations). In this work we have 
used only the longest of the exposures performed with the 
0° retarder plate angle. Science data and the associated 
calibrations have been retrieved through the ESO archive 
interfacell and reduced following the standard approach 
for de-biasing and fiat-fielding. The master fiat field has 
been obtained from the median-combination of sky images 
obtained during twilight with no polarimetric optics in. 
No observations of polarimetric standard stars were per- 
formed during the night, but analysis of the polarimetric 
data taken in the context of the polarimetric calibration 
plan shows that the instrument is stable (Fossati et al 
2007). We have finally obtained four reduced images cor- 
responding to the four retarder plate angles. As an exam- 
ple, we show in Fig. 1 the reduced image taken with the 0° 
retarder plate angle. We have checked for the presence of 
extended polarisation features which could be spatially as- 
sociated with the X-ray nebula detected around the pulsar 
by Chandra (Pavlov. 2001) but we have found none, thus 
confirming the quick analysis briefly reported in Mignani 
et al. (2003). 

2.2. Data Analysis 

We have then used our reduced images to compute the 
linear, phase-averaged polarisation degree of the Vela pul- 
sar. For a given source, the linear polarization degree Pl 
is calculated as: 

PL = {pA+p^y^' (1) 



http://archive.eso.org/archive/eso_data_products.html 



^ http://archive.eso.org/ 
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Fig. 1. VLT FORSl image of the Vela pulsar field taken 
with the R Bessel filter and polarised optics, which split 
the source image in ordinary and extra-ordinary beams. 
The pulsar counterpart {R = 23.9) is highlighted by the 
circle. 

where, following the notations defined in Landi 
Degl'Innocenti et al. (2007): 

Pq = y and Pv = j 

and I, Q, U are the Stokes parameters defined as in 
Shurcliff (1962). Applying the recipe described in Bagnulo 
et al.(2006), the Pq and Py parameters are computed ac- 
cording to the following relations: 




a=0 



Pu = -< 
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+ /V 



a=45 



a=67.5 



(2) 



(3) 



where f° and /® are the source counts measured in the 
ordinary and in the extraordinary mode, respectively for 
each of the four retarder plate angles (0°, 22.5°, 45° and 
67.5°). In our case, the source counts have been measured 
via standard aperture photometry using the photometric 
package of IRAF and using optimised apertures for both 
the source and the background. For all retarder plate an- 
gles we have obtained photometric accuracies of ~ 8 %, 
or better. To minimise the contamination from the three 
brighter stars close to the pulsar position (see Fig. 1) we 
have also measured the source counts through PSF pho- 
tometry. However, the results obtained with the two ap- 
proaches do not differ significantly. 

By substituting the values of f° and measured for 
the Vela pulsar for all of the four retarder plate angles in 



Eqs. (2) and (3), we have thus obtained Pq 0.02% ± 
3.7% and Pu = -9.4% ± 4.0%, where the errors on Pq 
and Pu are given by Eq. (3) of Fossati et al. (2007). These 
values have been then corrected as explained in the FORS 
User Manual and in Fossati et al. (2007) to account for 
the chromaticity of the retarder waveplate. These values 
are measured in the detector reference frame, which was 
oriented -1-10° North to East (see Sect. 2.1). In order to 
refer them to the celestial reference frame we have applied 
a counter-rotation with respect to the instrument position 
angle, using Eq. (9) of Landi Degl'Innocenti et al. (2007) 
with X = -10°. We thus end up with P^ = 3.6 % ± 3.7% 
and P{j = -8.7% ± 4.0% where P^ and P{j represent 
the linear polarization measured using the North Celestial 
Meridian as reference direction. 



2.3. Results 

By substituting Pq and P{j in Eq. (1), we have finally 
obtained Pl = 9.4 ±4%, where the associated error has 
been calculated using Eq. (6) of Fossati et al. (2007). This 
relatively large error is justified by the fact that, owing 
to the faintness of the target {R = 23.9), the uncertainty 
on the background subtraction dominates the photomet- 
ric errors on f° and and ultimately those on Pq and 
P{j. We have checked that our measurement is not af- 
fected by systematics effects like, e.g. instrumental or the 
background polarisation of the supernova remnant. The 
FORSl instrumental polarisation has been carefully char- 
acterized by Patat & Romaniello (2006) who found a sig- 
nificant spurious polarisation (at the level of w 1.5%) close 
to the edges of the CCD, while close to the field center is 
of the order 0.03 %, or less. Since our target has been po- 
sitioned at the center of the CCD, our measurement is 
not affected by the instrumental polarisation. We have fi- 
nally evaluated the background polarisation produced by 
the supernova remnant by repeating our measurements on 
several objects of different brightness close to the pulsar 
position, and we have found that its effect is probably 
less than 1 %, i.e. well within our global error budget. A 
better characterization of the background polarization in- 
troduced by the supernova remnants will be possible only 
when higher signal-to-noise data are available. While our 
measured polarisation degree is consistent with the one of 
Wagner & Seifert (2000), who reported Pl = 8.5%, our 
associated error turns out to be much larger with respect 
to their quoted 0.8 %. The source of the large discrepancy 
in the error estimate is not clear, although we speculate 
that Wagner & Seifert underestimated the error due to 
the neglect of the background subtraction contribution to 
the global error budget, since 0.8% is approximately the 
error that we would obtain considering the photon noise 
as the only source of measurement indetermination. 
From the values of P{j and Pq we have computed the lin- 
ear polarization position angle 9, defined as the angle from 
the North Celestial meridian to the major axis of the po- 
larisation ellipse (Landi Degl'Innocenti et al. 2007). For 
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Pq > 0, the correct expression to conrpute 6 is (Landi 
Degl'Irrnocenti et al. 2007): 



9 — — arctan 
2 




where 9 counted eastward from the North Celestial merid- 
ian. For our values of and Pq we obtain 9 = — 34°±11° 
or, equivalently, 9 = 146° ± 11°. The error on 9 has 
been computed using Eq. (7) of Fossati et al. (2007). We 
have checked for a possible coincidence between the pul- 
sar phase-averaged optical linear polarisation position an- 
gle 9 and the position angle 9x of the axis of symme- 
try of the X-ray arcs and jets detected around the pul- 
sar by Chandra (Pavlov et al. 2001). These structures are 
most commonly interpreted in terms of particle outflows 
from the pulsar equatorial wind and of coUimated outflows 
along the pulsar spin axis (Pavlov et al. 2001; Helfand et 
al. 2001). Alternatively, they are interpreted in terms of 
particle outflows along the pulsar's magnetosphere open 
field lines and of coUimated outflows along the pulsar mag- 
netic axis (Radhakrishnan & Deshpande 2001; Deshpande 
& Radhakrishnan 2006). For 9x different estimates exist, 
e.g. 310° ± 1.5° (Helfand et al. 2001), 307° ± 2° (Pavlov 
et al. 2001), and 310.63° ± 0.07° (Ng & Romani 2004). 
It has also been noted (e.g. Helfand et al. 2001) that 9x 
apparently coincides with the position angle 0^ of the pul- 
sar proper motion. For 6*^, based on HST observations 
Caraveo et al. (2001b) give 307° ± 2°, while Dodson et 
al. (2003) with the VLBI give 301° ± 1.8°. Interestingly, 
our measurement of 9 is not very different (within « 2a) 
from the available measurements of 9x and 0^, suggest- 
ing an alignment between the optical linear polarisation 
direction, the axis of symmetry of the X-ray arcs and jets, 
and the pulsar proper motion. Although the chance occur- 
rence probabilities of this alignment are not negligible, it 
is tantalising to speculate about it as a tracer of the con- 
nection between the pulsar's magnetospheric activity and 
its dynamical interactions with the surrounding medium. 
More precise measurements of the pulsar optical polar- 
isation, possibly supported by still to come polarisation 
measurements in X-rays, will hopefully provide a further 
observational handle for more robust theoretical specula- 
tions. 

3. Discussion 

We have then compared our measurement of the Vela pul- 
sar phase-averaged optical polarisation degree with the 
expectations of various pulsar magnetosphere models. We 
have used the polarisation code of Dyks et al. (2004, 
hereafter DHR04) which can calculate the linear polar- 
isation degree for the emission region of the outer gap 
model (Cheng et al. 1986; Romani & Yadigaroglu 1995) 
and its two-pole version (Dyks & Rudak 2003, hereafter 
DR03). The code carefully takes into account the macro- 
scopic properties of the emission region (3D spatial ex- 
tent, differential aberration, propagation delays, cumula- 



tion of radiation emitted from different parts of the emis- 
sion region at the same pulse phase) but it ignores much of 
micro physics. For instance, the intrinsic polarisation de- 
gree Pj?^"*'' of emitted radiation is assumed to be the same 
within the entire emission region (it is an input parameter 
of the code) and the polarisation angle is given by projec- 
tion of "bulk" electron acceleration on the observer's sky. 
The acceleration is only determined by the curvature of 
magnetic field lines and the corotation of pulsar magneto- 
sphere (no gyration component). Within the emission re- 
gion of a given model, the radiation power is independent 
on the altitude (the same number of photons is emitted 
per centimetre of electron's trajectory in the observer's 
frame). In the magnetic colatitude the emissivity is lim- 
ited to a narrow range of the magnetic field line footprint 
parameter. 

3.1. Comparison with the outer gap model 

For the outer gap model we assumed that the emission 
region extends from the null charge surface (Romani & 
Yadigaroglu 1995) up to rmax = 1-5 Pic with p^ax = 
0.95 Pic (Pic = cProt/ZTT is the fight cyfinder radius, 
where Pj-ot is the pulsar rotational period and c is the 
speed of light), we considered magnetic field lines with 
footprint parameter 0.75 ^ rove ^ 0.95, and we weighted 
the emission in magnetic colatitude by a Gaussian cen- 
tred at r^y^ = 0.85 (ct = 0.05). All the previous parame- 
ters are defined in DHR04. For the Vela pulsar emission 
geometry different estimates exist, e.g. a — 71°, C = 65° 
(Radhakrishnan & Deshpande 2001), where a is the dipole 
inclination and C is the viewing angle (see, e.g. Fig. 1 of 
DR03), C = 63.6° (Ng & Romani 2004) from two differ- 
ent geometrical models of the Vela pulsar X-ray nebula; 
a ~ 70°, C ^ 61° (DR03), from the two-pole model of 
the 7-ray profile; a ~ 55°, C ~ 49° (Johnston et al. 2001) 
from a fit to the position angle curve at 0.6 and 1.4 GHz; 
a - 137°, C ^ 143.5° (Johnston et al. 2005) from a fit 
at 1.4 GHz; ( = 53.8° (Krishnamohan & Downs 1983) 
from a fit of the position angle curve at 2.3 GHz and for 
a = 60°. In the following, we will assume a = 71° as a 
representative value. 

Fig. [2] shows the modelled Vela pulsar optical pulse 
profile and the phase-resolved polarisation properties for 
the representative case of a = 71° and ^ = 65° and for 
the same emission region parameters assumed above. As 
it is seen from the top panel, this set of parameters can 
reasonably well reproduce the observed optical pulse pro- 
file of the Vela pulsar (see, e.g. Gouiffes 1998), with two 
peaks separated by ~ 0.25 in phase. Thus, we are confi- 
dent that they provide a good approximation of the Vela 
pulsar's emission geometry. The two lower panels show 
the modelled polarisation position angle ip, measured from 
the projected pulsar rotational axis, and the depolarisa- 
tion factor DP = pmod/pmtr^ defined as the ratio between 
the modelled polarisation degree P^°'^ at the observer's 
location and the intrinsic polarisation degree of the emit- 
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Fig. 2. Phase-resolved polarisation properties of the Vela 
pulsar calculated for a simplified outer gap calculation 
(DHR04) with the emission region parameters described 
in the text and assuming a = 71° and C. = 65°. The pan- 
els show: a) the pulse profile; b) the modelled polarisation 
position angle ip; c) the depolarisation factor (see text for 
the definitions). The radio pulse is expected near phase 
zero. One and a half rotation period is shown. The light 
curve is over plotted in panels b and c for reference (grey 
line). 



ted radiation P^"'"'. From our model we have derived the 
phase-averaged value of the depolarisation factor, where 
the average is done over the Stokes parameters. We ob- 
tained DP = 0.68, a typical value for many sets of model 
parameters. From this value we have then computed the 
phase-averaged modelled polarisation degree P™°'^ for dif- 
ferent values of the intrinsic polarization P^"*"^- 
Fig. [3] presents the modelled phase-averaged polarisation 
degree P™°'^ of the Vela pulsar as a function of the view- 
ing angle C for two different values of the intrinsic optical 
polarisation p™*'' (see caption). As it is seen, for ( ~ 65° 
the modelled polarisation degree is much larger than the 
observed one if P™'"' = 75%^ while consistency would re- 
quire P™*'' to be as low as 13%. Fig. |3] also shows that 
for large a the modelled polarisation degree P™°d does 
not vary significantly for viewing angles ( ^ 40%. Thus 
we find, in agreement with Chen et al. (1996), that the 
phase-averaged polarisation Pl for the outer gap models 



(D 
■D 
C 
O 

a 

N 

o 

Q. 
T3 

to 



CO 

a 




Fig. 3. Modelled phase-averaged linear polarisation de- 
gree P™°'^ of the Vela pulsar as a function of the viewing 
angle ( for a simplified outer gap calculation (DHR04) and 
for a dipole inclination a = 71° (solid). As a reference, 
curves for a = 60° (dashed) and a = 80° (dotted) are also 
shown. The thick and the thin lines correspond to an in- 
trinsic optical polarisation P,^"*"' = 75% and P^*' = 13%, 
respectively. The range of viewing angles C derived with 
various methods (see text) is marked with the vertical 
shaded band. The measured phase-averaged polarisation 
degree Pl (this work) is shown as the horizontal grey band. 



is usually more sensitive to the intrinsic properties of the 
emitted radiation than to the global geometric parame- 
ters (a, (). While this makes it difficult to constrain a 
and C from the comparison of the modelled and observed 
Pl it also implies that our conclusions are less affected by 
the uncertain knowledge of the pulsar's emission geome- 
try. Thus, our simulation based on the standard outer gap 
model overpredicts the value of the phase-averaged polari- 
sation, unless the intrinsic polarisation is much lower than 
expected. This is most likely due to the fact that the code 
of DHR04 does not takes into account the convolution of 
the single electron emission beam with the structure of 
the emission region/magnetic field, which can be an im- 
portant source of depolarisation, cf. Fig. 4b and Fig. 6c 
in Chen et al. (1996) and eq. (42) in Epstein (1973). 
For large a, and for ( > 30° the two-pole caustic model 
(DR03), which includes the strong emission from below 
the null charge surface, predicts roughly two times smaller 
Pl than the standard outer gap model. However, the 
two-pole models (including the "extended slot gap" of 
Muslimov & Harding 2004) can not reproduce the ob- 
served optical pulse shape of the Vela pulsar. This may 
just indicate different emission regions for the optical and 
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Fig. 4. Phase-averaged modelled position angle ■0 as a 
function of the viewing angle ^ for the same range of dipole 
inclinations a as in Fig. [3l The range of viewing angles 
C, derived with various methods (see text) is marked with 
the vertical shaded band. 

for the 7-ray band, where the two-pole geometry performs 
quite well. 

We have then derived the phase-averaged modelled po- 
larisation position angle ^, where, as before, the average 
is done over the Stokes parameters. Fig. [4] presents the 
modelled phase-averaged polarisation position angle -0 as 
a function of the viewing angle C for different values of 
the dipole inclination a. Only one set of curves is drawn 
since '0 does not depend on the intrinsic polarization de- 
gree P™*"^. As it is seen, for C, > 40° the value of ij) de- 
pends very little on a and C- We note that the actual value 
of modelled position angle ^ depends on which emission 
mechanism is actually responsible for the optical emission 
(e.g. synchrotron radiation or curvature radiation). In the 
case of Fig. IH ip was computed assuming that the emit- 
ted waves, after appropriate microscopic averaging that is 
not modelled, are polarized parallel to the projected bulk 
acceleration of electrons. This direction of polarization is 
characteristic of curvature radiation (Jackson 1962) and 
it is parallel to projected trajectory of electrons' bulk mo- 
tion (with no gyration). Instead, if the radiation is polar- 
ized perpendicular to the direction assumed in our code, 
as is the case for the synchrotron radiation (Rybicki & 
Lightman 1979) one must correct ij) by subtracting 90°. 
So the value of the modelled position angle should be 
understood "modulo 90°". For instance, for ^ ~ 65° the 
modelled values of tp is either ~ 98° or 8° for polarization 
parallel or perpendicular to the projected bulk accelera- 
tion, respectively. This means that the modelled phase- 



averaged polarisation vector is roughly perpendicular or 
aligned to the pulsar rotation axis. For C on the other 
hemisphere {(,' = 180° — Q the derived position angle -0 
would just change sign i.e. -(/"(C) = see Table 1 in 

DHR04. 

A striking result visible in Fig. |4]is the proximity of "0 
to ~ 90° within the large range of a and C- To understand 
it we have recalculated the model of Fig. [5] with various 
kinematic effects temporarily removed from our code. We 
have run the code successively without the sweepback of 
magnetic field lines, without centrifugal acceleration and 
without aberration. This exercise has shown that it is the 
last effect (the aberration of photon emission direction as- 
sociated with transition from the corotating frame to the 
observer's frame) which is the most important in fixing 
■0 near ~ 90°. This contribution of corotational velocity 
to the emission direction in the observer's frame is large 
because the optical emission of the outer gap model orig- 
inates from altitudes that are a considerable fraction of 
the light cylinder radius. 

It is now interesting to compare the modelled polari- 
sation position angle ip with the observed one 9. However, 
we have to caution that a straight comparison between 
the two angles is not possible since the former is measured 
from the projected pulsar rotational axis, while the latter 
is measured from North. The modelled and observed po- 
larisation position angles are thus related by the equation 
= 9 — ipo, where tpo is the position angle of the pulsar's 
projected rotational axis defined in the usual way, i.e. mea- 
sured eastward from North. Thus, equating 9 = 146° ±11° 
(see Sect. 2.3) and ip allows for a direct measurement of V'o- 
For -0 = 98° ± 5° (polarization parallel to the bulk acceler- 
ation) we obtain -00 = 48° ± 12°. Instead, for = 8° ± 5° 
(polarization perpendicular to bulk acceleration) we ob- 
tain 00 = 138° ± 12°. The first case implies that the pul- 
sar projected rotation axis is almost perpendicular to the 
axis of symmetry of the X-ray jets and arcs, while the 
second case implies that it is almost aligned to them (see 
Sect. 2.3). If we assume that this geometry is the favoured 
one (e.g. Helfand et al. 2001; Pavlov et al. 2001), then 
^0 = 138° ± 12°, = 8° ± 5° and the optical radiation ac- 
tually is polarised perpendicular to the bulk acceleration 
assumed in our code. This is consistent with a synchrotron 
radiation origin of the optical emission. The other way 
around, if we assume that the optical emission is due to 
synchrotron radiation, our observed polarisation position 
angle would imply that the pulsar rotation axis is aligned 
with the axis of symmetry of the X-rays jets and arcs. 

Thus, our results would imply that there is a substan- 
tial alignment between the optical polarisation position 
angle, the pulsar rotation axis and proper motion vector, 
and the axis of symmetry of the X-ray nebula. Evidence 
for a similar alignment has been also found from phase- 
resolved radio polarisation measurements of a number of 
radio pulsars, including Vela (Johnston et al. 2005). 
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3.2. Comparison with other models 

We have qualitatively evaluated other models by compar- 
ing the predicted vs. observed polarisation properties pub- 
lished for the Crab pulsar which has a very similar geome- 
try to the Vela pulsar. The viewing angle C ~ 60° for both 
Crab and Vela (Ng & Romani 2004) and the dipole incli- 
nation a is believed to be large in both cases (e.g. Rankin 
1990). Furthermore, as we have shown above, the polari- 
sation properties are found to depend rather weakly on a 
and 

One of the more sophisticated models is the striped wind 
model (Petri & Kirk 2005). However, this model does 
not account for magnetic field irregularities in the emis- 
sion/reconnection region and this is probably why the 
phase-averaged polarisation degree calculated for the Crab 
pulsar is much larger (19-31%, Petri & Kirk 2005) than 
the observed on^. A version of the outer gap model by 
Chen et al. (1996) is also notably advanced, as it includes 
convolution of the single electron emission beam with spa- 
tial distribution, electron energy distribution and pitch 
angle distribution. The neglected magnetic field irregular- 
ities are probably less important in this case because the 
emission originates from a region well inside Ric. However, 
the assumed simplified geometry of the emission region 
(2D, no radial extent) seems not to be adequate as it does 
not allow to reproduce the observed pulse profile of the 
Crab pulsar. Accordingly, this model does not account for 
the depolarisation effect produced by the cumulation of 
radiation emitted from different altitudes. Using Fig. 6 
in Chen et al. (1996) one can calculate a partial phase- 
averaged polarisation degree that is limited to the phase 
range within ±44° around the maximum of main pulse. 
The result (9.4%) is two times larger than the value ob- 
served for the Crab pulsar within the same phase range 
(5.1% or 4.2% if the DC emission is excluded). 

4. Conclusions 

We have studied the Vela pulsar phase-averaged optical 
polarisation properties by reanalysing the data set used 
by Wagner & Seifert (2000). The fraction of linear, phase- 
averaged polarisation (Fl = 9.4% ± 4%) is qualitatively 
similar to what they found but with a larger error that we 
attribute to a more realistic error analysis, thus underscor- 
ing the need for observations with higher signal-to-noise 
ratio. In addition, we have measured for the first time the 
optical linear polarisation position angle {9 — 146° ± 11°) 
finding that it is apparently aligned with the axis of sym- 
metry of the X-ray arcs and jets and with the pulsar's 
proper motion vector. We have compared our results with 
the expectations of the outer gap as well as of various pul- 
sar magnetosphere models. Interestingly, for the outer gap 
model we have found that for synchrotron optical emis- 

^ In the case of the Crab pulsar, the observed phase-averaged 
polarisation degree is equal to 6.2%. After the "DC" compo- 
nent is subtracted (Kellner 2002; Kanbach et al. 2005) the 
phase averaged Pl becomes 3.7%. 



sion the observed linear polarisation position angle 9 im- 
plies that also the pulsar rotational axis is aligned with 
the axis of symmetry of the X-ray arcs and jets. For all 
models we found that the computed phase-averaged polar- 
isation degree is typically much larger than the observed 
one. The most probable reason for this is that these models 
are too idealised and do not account for some depolarisa- 
tion effects. Measurements of phase-averaged polarisation 
of pulsars can thus help to identify the weaknesses of vari- 
ous models and stimulate the development of the complex 
numerical codes required to calculate the expected polar- 
isation properties. More optical polarisation observations 
of pulsars, now in progress with the HST, will help to pro- 
vide the very much needed observational grounds to settle 
theoretical models. 
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